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ABSTRACT 


This  thesis  contains  a  detailed  experimental  and  numerical  investigation  of 
second-generation,  controlled-diffusion  compressor-stator  blades  at  an  off-design 
inlet-flow  angle  of  39.5°.  Investigation  of  the  blades  took  place  in  a  low-speed 
cascade  wind  tunnel  using  various  experimental  procedures.  The  objective  of  the 
wind  turmel  study  was  to  characterize  the  flow  field  in  and  around  the  blades  at 
the  off-design  angle,  and  to  investigate  flow  separation  near  the  mid-chord  for  a 
high  Reynolds  number  of  640,000.  It  was  known  from  previous  studies  that 
boundary  layer  thickness  on  the  end  walls  were  of  different  thicknesses.  Thus, 
prior  to  taking  data,  an  adjustment  to  the  end  wall  boimdary  layer  thickness  was 
attempted  by  insertion  of  an  aluminum  trip  strip  far  upstream  of  the  blades. 
Rake  probe  survey's  were  performed  upstream  and  downstream  of  the  blades  in 
order  to  obtain  spanwise  upstream  and  downstream  total  pressure  profiles. 
Surface  flow  visualization  was  performed  on  the  blades  using  a  titanium  dioxide 
and  kerosene  mixture.  Blade  surface  pressure  measurements  were  obtained 
using  a  40-hole  instrumented  blade  from  which  coefficients  of  pressure  were 
calculated.  A  standard  optics,  two<omponent  laser-Doppler  velocimeter  was 
used  to  characterize  the  flow  field  upstream,  in  the  boundary  layer  on  the  suction 
side  of  the  blades,  and  in  the  wake  region.  A  numerical  investigation  was 
conducted  using  the  rotor  viscous  code  3-D  developed  by  Dr.  Roderick  Chima  of 
NASA  Lewis  Research  Center. 

Overall,  good  agreement  between  flow  visualization,  blade  pressure 
measurements,  laser  measurements,  and  numerical  modeling  was  obtained. 
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I.  INTRODUCTION 


A.  BACKGROUND 

For  many  years,  jet  engines  have  been  limited  by  compressor  stall  and  off- 
design  performance  behavior.  Compressor  stall  can  lead  to  degradation  of 
engine  performance  and  possibly  the  loss  of  the  engine.  Research  and 
development  over  the  past  few  decades  has  worked  towards  increasing  engine 
blade  performance.  With  new  advances  in  compressor  blade  design  technology, 
such  as  computational  fluid  dynamics  (CFD)  analysis,  the  goal  is  to  improve 
engine  performance  by  allowing  higher  blade  loading  while  still  maintaining 
stall  margin  and  efficiency.  For  these  reasons,  a  new  generation  of  controlled- 
diffusion  (CD)  blading  was  developed. 

ControUed-diffusion  (CD)  blading  is  shaped  such  that  higher  angles  of 
incidence  may  be  achieved  before  boundary  layer  separation  occurs,  thus 
increasing  the  blade  loading.  This  is  done  by  designing  the  blade  to  control  the 
diffusion  on  the  suction  side  so  as  to  avoid  boundary  layer  separation.  Higher 
blade  loading  will  allow  more  turning  of  the  air  flow  for  a  given  number  of 
blades  (or  solidity),  or  the  same  turning  with  fewer  blades  (lower  solidity). 
Therefore,  fewer  blades  will  be  required  to  produce  the  same  compression  ratio 
which  will  result  in  a  lower  engine  weight  and  better  performance. 

The  present  study  conducted  at  the  Naval  Postgraduate  School  (NFS)  low- 
speed  cascade  wind  tunnel  (LSCWT)  involved  the  CD  compressor  stator  blades 
67B,  designed  by  Thomas  F.  Gelder  of  NASA  Lewis  Research  Center  [Ref.  1]. 
stator  67B,  together  with  rotor  67,  comprise  compressor  stage  67B.  The  67B  stator 
blades  were  second-generation  CD  blades  designed  as  an  improvement  over  the 
former  67A  first-generation  blades,  designed  by  Nelson  Sanger  [Ref.  2].  Prior  to 
the  study,  ten  midspan  stator  67B  compressor  blades  were  machined  from 
aluminum  and  installed  in  the  LSCWT. 
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Previous  studies  were  done  on  the  blades  at  a  design  inlet  flow  angle  of 
36.3°  by  Hansen  [Ref.  3],  and  at  an  off-design  inlet  flow  angle  of  38.0°  by 
Schnorenberg  [Ref.  4]. 

B.  PURPOSE 

The  objective  of  the  present  study  was  to  characterize  the  flow  pattern 
upstream,  in  the  passages  between  the  blades,  in  the  boundary  layer  of  the 
blades,  and  in  the  wake  region  at  an  inlet  flow  angle  of  39.5°  at  a  Reynolds 
number  of  640,000.  Also,  a  detailed  investigation  was  to  be  performed  on  the 
separation  region  which  occurred  near  mid-chord.  Various  methods  were  to  be 
used,  including  flow  visualization,  rake  probe  surveys,  blade  surface  pressure 
measurements,  and  laser-Doppler  velocimetry  (LDV)  measurements.  A 
adjustment  was  made  to  the  wind  tunnel  wall  boundary  layer  prior  to  taking 
data  to  obtain  improved  spanwise  symmetry  over  the  blade.  A  numerical 
investigation  using  3-D  CFD  was  also  conducted  and  the  predictions  were 
compared  with  the  experimental  results. 
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11. 


TEST  FACILITY  AND  SETUP 


A.  LOW-SPEED  CASCADE  WIND  TUNNEL 

The  low-speed  cascade  wind  tunnel  is  located  at  the  NFS  Turbopropulsion 
Laboratory  facilities.  The  wind  tunnel  is  powered  by  a  750  hp  electric  motor 
driving  a  turbo-vane  blower,  and  is  capable  of  producing  a  sustained  maximum 
free-stream  Mach  number  of  .4  in  the  test  section.  Figure  1  shows  a  schematic  of 
the  LSCWT  in  the  Low  Speed  Turbomachinery  Building  (Bldg.  213)  with  the 
associated  plenum  chamber,  drive  system,  and  inlet  and  exhaust  ducting. 
Hansen  [Ref.  3]  gave  a  detailed  description  of  the  test  facility  and  test  section. 
Tunnel  flow  conditions  were  documented  for  uniformity  and  periodicity  in  the 
cascade  test  section  using  20  Stator  67A  blades  at  approximately  40.0°  (design), 
43.0°  and  46.0°  inlet  flow  angle  by  Elazar  [Ref.  5]. 


Figiue  1.  NPS  Cascade  Wind  Tunnel  Facility. 
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B.  TEST  SECTION 

Figure  2  shows  the  layout  of  the  LSCWT  and  test  section  with  dimensions. 
Prior  to  the  present  study,  10  stator  67B  CD  blades  were  installed  in  the  test 
section  and  tested  at  the  design  inlet  flow  angle  of  36.3°  by  Hansen  [Ref.  3],  and 
at  38°  by  Schnorenberg  [Ref.  4].  Reference  4  contains  a  description  of  the 
procedure  to  adjust  the  inlet  flow  angle.  As  can  be  seen  in  Figure  2,  air  was 
forced  up  through  the  60  inlet  guide  vanes  were  it  was  turned  towards  the  test 
section.  The  flow  then  entered  the  test  section  and  was  once  again  turned 
vertically,  and  finally  exited  to  ambient  pressure. 


Figure  2.  Detailed  Schematic  of  Test  Section. 
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Each  blade  had  a  chord  length  of  127.25  mm  (5.01  in),  and  a  span  of  254 
mm  (10.0  in).  The  blades  were  separated  in  the  pitch-wise  direction  by  152.4  mm 
(6.0  in).  Figure  3  below  shows  a  detailed  profile  of  the  67B  stator  blade.  Test 
section  data  are  summarized  in  Table  1. 


Tunnel  Span 

254  mm  (10  in.) 

Blade  Type 

Stator  67B  Controlled-Diffusion 

Blade  material 

Aluminum 

Number  of  Blades 

10 

Blade  spacing 

152.4  mm  (6  in.) 

Chord 

127.14  mm  (5.01  in.) 

SoUdity 

0.834 

Thickness/Chord 

0.05 

Setting  Angles 

16.3°  ±0.1° 

Table  1.  Test  Section  Data. 
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Two  partially  instrumental  blades  containing  8  pressure  taps  each,  were 
installed  at  locations  2  and  8  (Fig.  2),  while  a  third  fully-instrumented  blade 
containing  42  pressure  taps  was  installed  at  location  6.  Blades  3  and  4  were 
treated  with  a  black  anodized  coating  to  minimize  laser  light  scatter  for  LDV 
measurements  and  were  also  used  for  flow  visualization. 

Eight  survey  stations  based  on  axial  chord  length  that  were  used  for  the 
study.  Figure  4  shows  the  distances  from  the  leading  edge  as  a  fraction  of  the 
blade  axial  chord.  The  station  numbering  system  was  consistent  with  earlier 
studies. 


Figure  4.  Survey  Stations  and  Numbering  in  Terms  of  Axial  Chord. 
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1.  Wall  Boundary  Layer  Adjustment 

It  was  shown  by  Webber  [Ref.  6]  using  upstream  rake  probe  surveys  that 
the  boundary  layer  thickness  on  the  north  wall  of  the  wind  tunnel  was  thicker 
than  on  the  south  wall.  At  high  Reynolds  number,  Schnorenberg  [Ref.  4] 
recorded  flow  separation  with  3-D  vortices  that  were  not  symmetrical  with 
respect  to  blade  chord  length.  Furthermore,  the  core  flow  was  displaced  towards 
the  south  wall  by  the  north  wall  boundary  layer,  making  the  flow  non- 
symmetrical  about  the  midspan  of  the  blade. 

The  objective  of  the  wall  boundary  layer  adjustment  was  to  improve  the 
flow  symmetry  about  the  midspan  of  the  blade.  A  1.5875  millimeter  (l/16m.) 
thick,  by  15.875  millimeter  (5/8in.)  wide,  by  1524  millimeter  (60in.)  aluminum 
trip  strip  was  inserted  into  the  flow  on  the  south  wall  just  upstream  of  the  60 
inlet  guide  vanes.  Figure  5  shows  the  schematic  of  the  aluminum  trip  strip 
inserted  into  a  holding  frame.  The  strip  caused  a  thicker  boundary  layer  to  form 
on  the  south  wall,  thus  displacing  the  flow  back  to  be  more  symmetrical  about 
the  midspan.  By  doing  so,  the  endwall  corner  vorticies  were  also  found  to  be  of 
the  same  size,  and  at  the  same  distance  from  the  leading  edge. 


South  side  of  tunnel 


! - 

■ 

238.125  mm 

(9375  in.) 

!■■■ 

1524  mm  (60  in.) 


Figure  5.  Aluminum  Boundary  Layer  Trip  Strip  in  Wooden  Frame. 
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C.  INSTRUMENTATION  SETUP 


1.  Flow  Visualization 

Blade  surface  flow  visualization  was  performed  using  a  titanium  dioxide 
(Ti02)  and  kerosene  mixture.  The  flow  patterns  on  the  surface  of  the  blades  were 
recorded  using  an  8mm  video  camera.  Black  and  white  photographs  were  taken 
after  the  solution  dried  on  the  blades  to  show  the  final  results.  The  steps  for 
mixing  the  solution  are  in  given  Appendix  A. 

2.  Blade  Siuface  Pressure  Measiurements 

Surface  pressure  measurements  were  recorded  using  a  48  channel 
pneumatic  Scanivalve  rotary  system  controlled  by  a  HP-9000  computer.  The 
software  to  control  the  system  was  fully  documented  by  Classick  [Ref.7],  and 
later  modified  by  Armstrong  [Ref. 8].  Scanivalve  ports  and  channel  assignments 
are  shown  in  Table  Cl  in  Appendix  B. 

3.  Rake  Probe  Measurements 

A  20-hole  rake  probe  was  used  to  acquire  pitch-wise  surveys  of  the 
spanwise  distribution  of  coefficient  of  pressure  Cp  upstream  and  downstream  of 
the  blades.  The  rake  probe  consisted  of  17  total  pressure  ports,  1  static  pressure 
port,  and  2  yaw  ports.  Figure  6  shows  a  schematic  of  the  probe  along  with  an 
upstream  survey  diagram  showing  position  and  traverse  direction.  Here  again, 
data  were  recorded  using  a  48  chcinnel  pneumatic  Scanivalve  rotary  system 
controlled  by  a  HP-9000  computer.  Scanivalve  ports  and  channel  assignments  are 
shown  in  Table  C2  in  Appendix  B. 
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Figure  6.  Upstream  Survey  Schematic  and  Rake  Probe. 

4.  LDV  Measurements 

LDV  measurements  were  obtained  using  a  four-beam,  two-color  TSI 
Model  9100-7  system.  A  description  of  the  setup,  including  optics,  atomizer 
seeding,  and  laser  type  are  thoroughly  described  by  Elazar  [Ref.  5].  Data 
acquisition  and  the  traverse  mechanism  were  controlled  by  a  personal  computer 
(PC)  using  TSI  Flow  Information  Display  (FIND)  Software  (version  4.0). 
Software  inputs  used  for  the  present  study  are  given  in  Appendbc  C.  The  data 
acquisition  and  traverse  mechanism  are  described  in  detail  by  Murray  [Ref.  9]. 
The  LDV  laser  and  traverse  table  are  shown  in  Figure  7. 
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Figure  7.  LDV,  Traverse 


and  Cascade  Test  Section. 


III. 


EXPERIMENTAL  PROCEDURE 


A.  WALL  BOUNDARY  LAYER  ADJUSTMENT 

Various  methods  such  as  meshed  screens  (with  different  types  of 
fineness),  honeycomb  screens,  and  flow  straighteners  have  been  used  to  adjust 
inlet  flows  and  to  make  them  more  uniform.  The  procedure  tried  here  was  to 
insert  various  widths  of  aluminum  strips  into  the  flow  field  on  the  south  side  of 
the  wind  tunnel  wall,  just  upstream  of  the  60  inlet  guide  vanes  so  as  to  make  the 
flow  on  the  blade  surface  more  S5nnmetrical  in  the  spanwise  direction.  After 
installation  of  ihe  aluminum  boundary  layer  trip  strip,  a  titanitim  dioxide  and 
kerosene  solution  was  brushed  onto  blades  3,  4,  and  5.  The  wind  tunnel  was 
started  and  allowed  to  stabilize  at  a  test  section  speed  of  Mach  .22  (s  70.0  m/s). 
Surface  flow  visualization  was  used  to  record  the  flow  transient  and  steady-state 
flow  patterns  on  the  blades.  This  procedure  was  done  repeatedly  for  various 
widths  of  aluminum  strips.  It  was  found  that  a  15.875  mm  protrusion  of  an 
aluminum  strip  into  the  airflow  gave  symmetrical  flow  about  the  blades,  which 
appeared  to  be  symmetrical  in  the  spanwise  direction. 

B.  WIND  TUNNEL  CALIBRATION 

During  the  calibration  runs,  the  wind  tunnel  was  allowed  to  reach  a 
plenum  temperature  equilibrium  for  each  speed.  The  tunnel  was  run  at  7 
different  speed  settings  (plenum  pressure),  78.74  mm  (3.1  in.),  114.3  mm  (4.5  in.), 
147.32  mm  (5.8  in.),  203.2  mm  (8.0  in.),  254  mm  (10.0  in.),  304  mm  (12.0  in.),  and 
355.6  mm  (14.0  in.)  of  H2O.  Plenum  pressure,  plenum  temperature,  and  ambient 
pressure  were  recorded.  Using  the  LDV,  horizontal  (U)  and  vertical  (V)  velocities 
components  were  recorded  for  each  speed.  The  data  were  used  in  a  FORTRAN 
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program ,  CALIB1.FOR,  which  fit  the  tunnel  characteristics  using  a  least-squares 
method  to  determine  the  pressure  ratio  as  a  function  of  Mach  number  [Ref.  3]. 

C.  FLOW  VISUALIZATION 

Surface  flow  visualization  was  performed  using  a  titanium  dioxide  (Ti02) 
and  kerosene  mixture.  Mixing  procedures  are  given  in  Appendbc  A.  With  the 
acrylic  window  removed,  the  mixture  was  applied  evenly  on  blades  3,  4,  and  5 
with  a  fine-hair  paint  brush.  The  acrylic  window  was  immediately  reinstalled 
and  the  wind  tunnel  was  started.  The  tunnel  was  bought  up  to  a  speed  of  Mach 
•22  in  the  test  section.  An  8mm  VHS  video  camera  mounted  on  a  tripod  was  used 
to  record  the  transient  and  final  flow  field  patterns  on  the  blades.  All  test  section 
flow  conditions  corresponded  to  a  Reynolds  nimiber  of  640,000. 

D.  LDV  MEASUREMENTS 

1.  Probe  Volume  Alignment 

Prior  to  each  survey,  the  probe  volume  formed  by  the  intersecting  laser 
beams  was  aligned  with  an  aluminum  alignment  tool.  Details  on  alignment 
procedures  are  described  in  reference  4.  All  surveys  were  done  at  midspan  of  the 
blades. 

2.  LDV  Surveys 

The  LDV  was  aligned  and  leveled  such  that  the  X,  Y,  Z,  traverse  motions 
would  move  the  measurement  volume  horizontal  (blade-to-blade),  vertical 
(normal  to  the  leading  edge  locus),  and  parallel  (spanwise  to  the  blade),  for 
surveys  taken  at  stations  1  through  13.  For  station  3,  the  laser  was  pitched  up  5°, 
while  at  station  11,  the  laser  was  pitched  down  5°,  in  order  to  avoid  interference 
of  the  blue  beam  with  the  leading  and  trailing  edge  respectively.  For  all 


12 


boundary  layer  surveys,  the  laser  was  yawed  4°  to  the  left  in  order  to  avoid 
interference  of  the  green  beam  with  the  end  tip  of  the  blade.  A  total  of  12  LDV 
survey's  were  done  at  an  off-design  inlet  flow  angle  of  39.5°  for  a  Reynolds 
number  of  640,000.  Boundary  layer  surveys  were  completed  on  blade  #3  with 
station  8  repeatability  measurements  on  blade  #4.  Inlet  flow  surveys  were  done 
at  stations  1  and  3,  while  wake  survey's  were  done  at  stations  11  and  13.  With  the 
measurements  taken  in  coincidence  mode,  a  total  of  1000  data  points  were 
collected  for  each  sample  point.  Axial  (vertical)  velocities  U,  were  recorded  using 
the  514.5  nm  green  beams,  while  tangential  (horizontal)  velocities  V,  were 
recorded  using  the  488  nm  blue  beams.  Fringe  spacing  based  on  half-angle 
calculations  gave  4.7569  microns  for  the  green  beam ,  and  4.5119  microns  for  the 
blue  beam.  A  5  Mhz  frequency  shift  was  used  to  detect  flow  reversal. 

For  each  survey,  plenum  total  pressure  (Pto)  ,  plenum  total  temperature 
(Tto),  and  ambient  pressure  (Pamb)  were  recorded.  Program  C  ALIBI  .FOR 
[Hansen,  Ref.  3],  used  Pt©,  Tto,  and  Pamb  to  calculate  the  tunnel  inlet  reference 
velocity  (Vref)  for  each  survey.  Vref  is  then  used  to  non-dimensionalize  the  total 
velocity  (W),  axial  velocity  (U),  tangential  velocity  (V),  and,  U  and  V  turbulence 
intensities,  so  that  individual  surveys  could  be  compared.  Appendix  D  lists  the 
non-dimensionalized  data  for  each  survey. 
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IV. 


RESULTS  AND  DISCUSSION 


A.  FLOW  VISUALIZATION 

Results  show  that  after  the  boundary  layer  trip  strip  was  inserted  into  the 
upstream  flow  field  on  the  south  wall  of  the  tunnel,  the  inlet  flow  was  more 
symmetrical  about  the  midspan  of  the  blades.  Figures  8  and  9  show,  for 
comparison,  Schnorenberg's  [Ref.  4]  results  and  the  present  results,  respectively. 
For  the  present  study,  two  counter-rotating  vortices  appeared  at  approximately 
•78  Cac.  .  These  vortices  were  the  result  of  corner  vortices  which  formed  due  to 
the  interaction  of  the  endwall  boundary  layer  with  the  blade  leading  edge.  Both 
vortices  were  s)mrimetric  about  the  spanwise  direction,  and  located  at  the  same 
chordwise  position. 

Measurements  taken  from  the  photographs  showed  that  separation  of  the 
flow  occurred  at  .5  Cac  at  the  nudspan  of  the  blades.  The  separation  line  was  not 
straight  along  the  span  of  the  blade  because  of  the  three-dimensional  endwall 
effects.  Actual  separation  of  the  boxmdary  layer  was  most  probably  farther 
downstream  because  of  the  gravitational  effects  of  the  titanium  dioxide  and 
kerosene  mixture. 
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Figure  8.  Previous  Result  without  Boundary  Layer  Trip  Strip. 


Figure  9.  Result  with  Boundary  Layer  Trip  Strip. 
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Flow  visualization  also  showed  excellent  periodicity  between  blades  3,  4, 
and  5.  Figure  10  is  a  photograph  taken  of  blades  3  and  4  to  show  how  periodic 
the  flow  was  from  blade  to  blade.  An  averaged  5%  difference  was  measured 
between  the  locations  of  the  vortices  and  separation  points  between  blades  3,  4, 
and  5. 


Figure  10.  Flow  Visualization  Periodicity  study. 
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B.  BLADE  SURFACE  PRESSURE  MEASUREMENTS 

Blade  surface  pressure  measurements  were  taken  on  blade  6  at  the  high 
Reynolds  number.  Figure  11  below  shows  the  results  of  the  pressure  distribution 
in  terms  of  the  coefficient  of  pressure,  Cp ,  vs  fraction  of  blade  chord  (^/ c). 


Figure  11.  Cp  vs  for  Beta  =  39.5°. 


For  the  suction  side  of  the  blade,  Cp  is  shown  to  decrease  to  -0.8  almost 
immediately  about  the  leading  edge.  The  diffusion  between  suction  point  2  and 
3  may  be  an  indication  of  a  leading  edge  separation  bubble.  Re-acceleration  of 
the  flow  continues  to  0.4  ^/c  where  a  minimum  Cp  of  -1.17  was  reached.  This 
point  corresponded  to  a  maximum  Mach  nximber  of  .345.  From  0.4  %/ c  the  Cp 
distribution  gradually  increased  linearly  to  0.64  ^/c  which  showed  no  sign  of 
flow  separation.  However,  between  0.64  and  0.69  ^/c  a  severe  adverse 
pressure  gradient  existed  causing  turbulent  flow  separation  on  the  blade.  The  Cp 
distribution  over  most  of  the  pressure  surface  was  constant,  except  at  the  trailing 
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edge.  The  jump  in  was  caused  by  a  reverse  flow  aft  of  the  blunt  trailing  edge 
of  the  blade. 

A  comparison  was  made  of  the  on-design  and  two  off-design  angles  of 
incidence.  On-design  blade  pressure  measurements  were  obtained  by  Hansen 
[Ref.  3]  at  36.3°,  while  Schnorenberg  [Ref.  4]  performed  measurements  at  38.0° 
inlet  flow  angle.  Figure  12  shows  the  comparison  of  the  results.  Overall,  they 
compared  well.  The  current  study,  at  39.5°,  showed  that  for  the  region  between 
0.0  and  0.4  ^/c  ,  a  higher  blade  loading  occurred  with  a  possible  separation 
bubble  at  the  leading  edge.  After  0.4  ^c  a  slightly  higher  diffusion  rate  was 
seen  with  strong  diffusion  between  0.64  ^/c  and  0.69  ^/c.  No  significant 
differences  were  noted  for  the  pressure  side  of  the  blade. 


Figure  12.  Comparison  of  Cp  vs.  Plots  for  P  =  36.3°,  38.0°,  and  39.5°. 
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C.  RAKE  PROBE  MEASUREMENTS 

An  upstream  inlet  survey  across  passage  3  was  made  using  a  20-hole 
pneumatic  rake  probe.  Total  and  static  pressure  measurements  were  obtained 
and  plotted  as  Cp  vs  tunnel  span,  as  shown  in  Figure  13  0).  An  averaged  inlet 
boundary  layer  thickness  was  foimd  to  be  between  50.8-69.85  mm  (2.0  -  2.75  in.) 
on  the  north  wall  of  the  tunnel,  and  about  76.2  mm  (3  in.)  on  the  south  wall.  Non- 
uniform  boundary  layers  were  attributed  to  wakes  from  the  inlet  guide  vanes  not 
being  fully  mixed  out,  particularly  in  the  endwall  regions.  The  Cp  distribution  in 
the  center  of  the  tunnel  was  constant  over  40%  of  the  span. 


Figiure  13.  Upstream  Spanwise  Cp  Distribution. 

(1)  Here,  Cp  is  defined  as  the  non-dimensional  pressure,  i.e.,  the  local  total  pressure  divided  by  the 
maximum  total  pressure  measured  by  the  rake  probe. 
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D.  LDV  MEASUREMENTS 


From  blade  surface  flow  visualization,  it  was  shown  that  air  flow  was 
mostly  symmetrical  at  the  midspan.  Separation  of  the  boundary  layer  occurred  at 
0.5  Cac  The  Cp  plot  obtained  from  the  fully  instrumented  blade  showed  that  a 
possible  separation  occurred  between  0.64  Cac  and  0.69  Cac .  LDV  meastuements 
at  the  midspan  section  were  made  for  comparison  with  previous  experimental 
data  and  to  obtain  a  more  detailed  picture  of  the  flow  field  characteristics  around 
the  blade. 


1.  Inlet  Surveys 

a.  Station  1 

Results  showed  nearly  uniform  velocity  ratio's  W/Vref ,  U/Vref  , 
and  V/Vref  as  shown  in  Figure  14.  The  wave-like  features  of  the  velocity  ratios 
were  caused  by  upstream  influence  of  the  blade  profiles  since  the  disturbance 
repeated  itself  every  blade  spacing.  Axial  and  tangential  velocities  gave  an 
average  inlet  flow  angle  of  39.5°  with  a  turbulence  intensity  of  2%  for  both  the  U 
and  V  velocity  components.  There  was  some  indication  of  the  unmixed  inlet 
guide  vane  wakes  in  the  turbulence  intensity  data  for  the  tangential  velocity  (Tv) 
component  which  showed  a  three-per-blade  spacing  ripple.  This  could  be  due  to 
wakes  that  had  coupled  together,  however,  this  effect  was  not  repeatably 
measured.  The  Reynolds-stress  correlation  coefficient  remained  at  a  constant 
value  of  0.1,  indicating  a  random  or  uncorrelated  flow. 
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It  was  found  by  LDV  measurements  that  inserting  the  trip  strip 
caused  the  inlet-flow  angle  (Pi)  to  changed  from  an  average  38.0°  to  39.5°.  Since 
the  inlet-  flow  angle  and  side-wall  angle  were  not  the  same  (1.5°  difference), 
there  was  most  likely  a  mild  secondary  flow  entering  the  test  section  as  shown  in 
Figure  15.  Repeatability  tests  were  performed  at  this  station  with  good  results  on 
the  mean  flow  quantities. 


Figure  15.  Possible  Secondary  Airflow  Pattern. 


b.  Station  3 

The  survey  results  at  station  3  are  shown  in  Figure  16.  Velocity 
ratio's  W/Vref,  and  U/Vref  decreased  toward  the  leading  edge  of  each  blade.  This 
was  a  result  of  potential  effect  of  the  blades  on  the  approaching  flow.  The  V 
velocity  component  was  actually  accelerated  around  the  leading  edge,  thus  the 
increase  in  V/Vref  was  seen  as  the  flow  approached  the  blade.  Turbulence 
intensity  for  both  U  and  V  velocity  components  stayed  at  2.0%  with  a  slight 
increase  to  2.5  %  on  the  suction  side  of  the  leading  edge;  while  the  correlation 
coefficient  remained  at  0.1. 
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Velocity  Ratio 


2.  Boundary  Layer  Surveys 

Boundary  layer  surveys  were  performed  at  stations  6,  7,  8,  and  9  on  the 
suction  side  of  blade  #3.  A  test  for  periodicity  was  performed  at  station  8  on 
blade  #4.  Boundary  layer  data  are  presented  in  terms  of  the  non-dimensional 
distance  (d/c)  perpendicular  to  the  blade  surface,  where  c  is  the  blade  chord 
length.  The  distance  between  each  data  point  taken  was  0.5  mm  perpendicular  to 
the  surface  of  the  blade  for  all  stations.  Comparison  with  previous  results  are 
discussed. 


a.  Station  6 

The  results  obtained  at  station  6  (0.25  Cac),  are  shown  in  Figure  17. 
Thirty  eight  data  points  were  taken  perpendicular  to  the  blade  surface.  Results 
showed  that  the  flow  was  turbulent  and  attached  to  the  blade.  Acceleration  of  the 
flow  to  1.3  times  the  inlet  reference  velocity  was  measured.  The  second  measured 
point  (at  a  d/c  of  0.008)  was  2%  higher  than  the  first  point,  which  indicated  that 
the  edge  of  the  boundary  layer  was  approximately  at  that  distance  from  the 
blade  surface.  As  shown  in  Figure  17,  W/Vref ,  U/Vref ,  and  V/Vref  gradually 
decrease  in  magnitude  as  distance  from  the  blade  increased  and  the  pressure  side 
of  the  passage  was  approached.  Turbulence  intensities  remained  relatively 
constant  at  2.0  %  over  most  of  the  survey.  The  correlation  coefficient  was  0.1. 
Previous  experimental  LDV  surveys  at  Pi=38.0°  (off-design)  at  640,000  Re, 
showed  similar  results  [Ref.  4]. 
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b. 


Station  7 


Station  7  (0.50  Cac),  was  surveyed  by  collecting  40  data  points 
perpendicular  to  the  blade  surface.  Results  are  shown  in  Figure  18.  The  initial 
increase  in  W/Vref ,  U/Vref ,  and  V/Vref  was  within  the  boundary  layer,  where 
the  thickness  was  measured  at  0.02  d/ c  (2.54  mm  out  from  the  blade  surface). 
The  maximum  total  velocity  ratio  at  the  edge  of  the  boundary  layer  was  1.1. 
Tmbulence  intensity  for  the  U  velocity  reached  a  maximum  of  18.0%  at  the 
surface  of  the  blade  and  decreased  to  2%  at  the  outer  edge  of  the  boundary  layer. 
Turbulence  intensity  for  the  V  velocity  component  remained  constant  at  2% 
throughout  the  survey.  The  correlation  coefficient  first  rose  from  a  value  of  0.0  to 
0.1  at  the  surface  of  the  blade  and  then  decrease  to  a  value  of  -0.41  at  the  outer 
edge  of  the  boundary  layer.  It  then  gradually  increased  to  a  value  of  0.0  at  the 
end  of  the  survey. 

The  survey  showed  that  the  boundary  layer  was  most  probably  still 
attached,  thus  meaning  that  the  separation  point  was  further  downstream  than 
this  station.  This  indicated  that  flow  visualization  results  were  contaminated  due 
to  gravitational  effects.  Comparing  the  velocity  ratios  obtained  here  with 
previous  results  showed  that,  for  the  38.0°  off-design  incidence  [Ref.  4],  the 
results  were  sunilar.  The  turbulence  intensities  were  also  similar. 
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Figure  18.  Boundary  Layer  LDV  Survey  at  Station  7. 


c. 


station  8 


Station  8,  (0.75  Cac)  ,  was  surveyed  by  collecting  43  data  points 
perpendicular  to  the  blade  surface.  Results  shown  in  Figure  19  indicated  that  air¬ 
flow  reversal  was  measured  within  the  boundary  layer  from  0.0  to  0.06  d/c.  The 
reverse  flow  velocity  magnitude  was  approximately  10%  of  the  reference 
velocity  (Vref).  The  V/Vref  ratio  maintained  a  constant  value  of  0.0  until  the  value 
of  0.06  d/ c  was  reached,  at  which  time  it  increased  gradually  to  a  value  of  0.3. 
This  indicated  that  the  V  component  velocity  vector  was  always  in  the  positive 
direction,  i.e.,  diverging  away  from  the  blade  surface. 

The  axial  turbulence  intensity  (Tu)  survey  followed  a  Gaussian 
distribution.  Initially  it  started  at  5%,  climbed  to  30%  at  0.085  d/c,  and  then 
dropped  back  down  to  5%.  The  maximum  turbulence  corresponded  to  the 
maximum  shear  gradient  in  the  axial  velocity  distribution.  Tangential  turbulence 
intensity  Tv,  only  ranged  from  5-9%.  The  correlation  coefficient  ranged  from  0.1 
to  -0.1  throughout  the  survey. 

A  periodicity  test  was  done  at  station  8  on  blade  #4  to  see  how  well 
data  matched  with  those  at  station  8  on  blade  #3.  Results  matched  exceptionally 
well  with  no  significant  differences,  as  can  been  seen  in  Figure  20.  Periodicity 
was  also  confirmed  with  the  flow  visualization  pictures,  (Fig.  10). 
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d. 


Station  9 


Station  9  (0.95  Cac)  ,  was  surveyed  by  collecting  58  data  points 
perpendicular  to  the  blade  surface.  Results  given  in  Figure  21,  showed  that  flow 
reversal  appeared  within  the  boundary  layer  from  the  blade  surface  to  a  distance 
of  0.1  d/ c.  This  can  be  seen  in  the  distribution  of  the  U/Vref  velocity  ratio.  From 
0.1  d/ c,  U/Vref  increased  to  a  maximum  value  of  1.0  at  the  outer  edge  of  the 
boundary  layer.  V/Vref  velocity  ratio  started  out  with  a  positive  value  of  0.1  and 
gradually  decreased  to  a  value  of  -0.1.  This  indicated  that  the  flow  direction  was 
first  positive  (to  the  right  of  the  vertical)  and  then  eventually  became  negative  ( 
to  the  left  of  vertical). 

Turbulence  intensity  Tu  ranged  from  5%  to  30%,  while  turbulence 
intensity  Tv  ranged  from  5%  to  20%.  The  correlation  coefficient  began  with  a 
value  of  0.0  and  gradually  climbed  to  0.3  at  0.125  d/c  and  then  gradually 
decreased  back  to  0.0. 

Comparing  the  current  study  to  the  off-design  case  at  Pi=38.0° 
[Ref.  4],  little  difference  was  foimd  in  the  results.  No  simularity  in  data  was 
evident  when  compared  to  the  on-design  case,  at  pi  =36.3°  [Ref.  3]. 
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3.  Wake  Surveys 

Wake  surveys  were  performed  at  stations  11  and  13  over  two  blade 
passages  as  shown  in  Figures  22  and  23  respectively.  Both  had  similar  results  and 
therefore  only  station  13  will  be  discussed. 

a.  Station  13 

Station  13  was  surveyed  across  two  passages  collecting  52  data 
points.  The  velocity  profiles  indicate  a  minimum  at  the  trailing  edge  of  each 
blade,  noted  by  the  deficit  in  the  axial  velocity  distribution.  In  the  freestream 
region,  the  velocity  ratios  indicated  a  slight  decrease  in  magnitude  on  the  suction 
side  of  the  blade,  dien  increased  slightly  as  they  approached  the  pressure  side  of 
the  blade.  The  average  exit  angle  was  calculated  to  be  9.5°  from  the  axial 
direction. 

The  axial  turbulence  intensity  showed  two  peaks  at  the  trailing 
edge  of  the  blades,  with  a  maximum  value  of  28%.  Wake  tangential  velocity  also 
had  two  peaks  with  a  maximum  value  of  20%.  The  correlation  coefficient  varied 
from  0.4  to  -0.4.  The  wake  thickness  was  approxinaately  3.3%  thicker  for  the 
current  study  when  compared  to  the  previous  study  [Ref.  4].  Comparing  the 
results  to  the  on-design  [Ref.  3]  case  showed  that  the  velocity  ratios  for  the  on- 
design  case  were  shifted  approximately  7.2  mm  to  the  left,  i.e.,  exit  flow  was 
turned  more  through  the  passage.  Average  exit  flow  angle  was  approximately 
1.0°  compared  to  9.5°  for  the  current  study.  Ftuthermore,  wake  thickness  was 
approximately  14%  less  than  that  of  the  current  study. 
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V. 


COMPUTATIONAL  FLUID  DYNAMIC  (CFD)  ANALYSIS 


A  PURPOSE 

The  purpose  of  the  numerical  analysis  was  to  obtain  a  solution  that  could 
be  validated  against  experimental  data  obtained  at  an  off-design  incidence  angle 
of  39.5°.  Specifically,  coefficient  of  pressures,  flow  reversal,  and  vortex  locations 
could  be  compared.  By  validating  the  CFD  solution  with  the  experimental  data, 
confidence  is  gained  in  the  use  of  the  code  to  arrive  at  designs  that  give 
improved  blade  performance. 

B.  NUMERICAL  PROCEDURE 

Computational  fluid  dynamic  (CFD)  analysis  was  performed  using  the 
Rotor  Viscous  Code  3-D  (RVC3D  -  version  920318)  developed  by  Dr.  Roderick 
Chima  of  NASA  Lewis  Research  Center  [Ref.  11].  RVC3D  is  a  computer  code  for 
analysis  of  three-dimensional  viscous  flow  in  turbomachinery.  The  code  solved 
the  thin-layer  Navier-Stokes  equations  with  an  explicit  finite-difference 
technique.  Txurbulence  effects  were  modeled  using  a  3-D  adaptation  of  the 
Baldwin-Lomax  turbulence  model.  The  equations  were  discretized  using 
second-order  finite-differences  and  were  solved  using  a  four-stage  Runge-Kutta 
scheme. 

A  3-D  grid  was  generated  around  a  single  blade  for  half  its  span.  A 
constant  Courant  ntunber  (CFL)  of  5.0  was  used  throughout  all  the  calculations. 
Grid  generation  and  RVC3D  inputs  for  numerical  analysis  can  be  found  in 
Appendix  E. 

The  Cp  distribution  from  the  inlet  rake  probe  survey  was  used  to  calculate 
the  inlet  boundary  layer  thickness  on  the  endwall.  The  test  section  Mach  number 
was  .22  at  an  inlet  air  angle  pi=39.5°.  Angle  of  incidence  was  changed  by 
changing  the  parameter  'prat'  (static  pressure  of  hub  exit  to  inlet  reference  total 
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pressure  (Phubexit  /Po)).  Prat  was  calculated  from  rake  probe  measurements  to  be 
0.9729  for  a  inlet  angle  of  39.5°. 

The  parameter  "jedge"  is  defined  as  the  last  j-index  (away  from  the  airfoil) 
searched  for  the  turbulent  length  scale.  For  the  Balwin-Lomax  turbulent  model, 
"jedge"  should  be  a  grid  line  slightly  bigger  than  the  largest  expected  blade 
boundary  layer.  Initially  this  was  set  to  the  maximum  of  49,  which  was  much 
bigger  than  the  boundary  layer  thickness  indicated  by  the  experimental  data. 

The  "kedge"  was  the  maximum  grid  surface  (away  from  the  endwall)  to 
which  the  flowfield  was  searched  for  the  turbulent  length  scale.  The  "jedge"  and 
"kedge"  parameters  were  relaxed  from  an  inital  value  of  49  and  70  to  a  value  of 
30  and  50  respectively.  One  other  parameter  varied  was  "cmutm",  the  value  of 
(eddy  viscosity)/  (laminar  viscosity)  at  which  transition  from  laminar  to 
turbulent  is  assumed  to  occur.  Typically  a  value  of  14  is  normally  used  for 
natural  transition  and  a  value  of  0.0  is  used  to  simulate  fully  turbulent  flow.  The 
number  was  reduced  to  10  because  the  flow  was  turbulent  from  close  to  the 
leading  edge.  Appendix  F  contains  the  output  for  the  inlet  and  exit  conditions 
that  span  in  the  K  direction. 

C.  GRID  GENERATION 

A  two-dimensional  ^id  was  first  computed  using  a  modified  version  of 
the  FORTRAN  code  GRAPE  (Grids  About  Airfoils  using  Poisson's  Equations). 
Reference  Hansen  [Ref.  3]  for  code  inputs.  The  grid  size  was  340  x  49.  Grid 
coordinates  were  generated  based  on  manufacturing  dimensions.  Next,  a  three- 
dimensional  grid  was  built  using  a  FORTRAN  program  called  STACK,  which 
took  the  two  dimensional  C-type  grid  and  extended  it  outward  in  the  z-direction 
for  70  grid  points.  The  final  grid  ended  up  being  a  340  x  49  x  70,  which  consisted 
of  approximately  1.2  million  grid  points.  Figure  24  shows  the  final  grid. 
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D.  RESULTS  AND  DISCUSSION 


1.  Density  Residual  History 

Figure  25  shows  the  density  residuals  up  to  seven  thousand  iterations.  The 
residuals  started  out  at  a  approximately  5.0  x  10'^  and  decreased  by  three  orders 
of  magmtude  in  7000  iterations.  It  took  approximately  30  seconds  per  iteration 
on  the  CRAY  J90  at  NFS. 


Density  Residuals  vs  Number  of  Iterations 


Number  of  Iterations 
Figure  25.  Convergence  History. 

2.  FAST  Flow  Analysis 

CFD  flow  visualization  results  showed  good  correlation  with 
experimental  studies.  Figure  26  shows  particle  traces  that  indicated  a  vortex  at 
approximately  the  same  location  as  that  of  the  experimental  flow  visualization. 
The  red  particle  traces  represent  flow  that  is  forced  from  the  endwall  towards  the 
midspan  of  the  blade.  As  can  be  seen  with  the  yellow  lines,  some  reverse  flow 
over  most  of  the  blade  and  vortex  formation  near  the  endwall  corner  occured. 
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occured.  CFD  analysis  did  not  pick  up  any  reverse  flow  directly  at  the  midspan, 
nor  any  separation  bubble  near  the  leading  edge. 


Figure  26.  Particle  Traces  of  the  Flow  Field  over  the  Suction  Surface. 
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3. 


Coefficient  Of  Pressure  Distribution 


Results  for  CFD  vs.  Experimental  Cp  distributions  are  shown  in  Fig. 
27.  There  was  an  immediate  decrease  in  pressure  around  the  leading  edge 
(suction  side)  of  approximately  -0.7  which  sharply  increased  to  -0.3  as  predicted 
by  the  code.  The  solution  was  closely  inspected  with  FAST,  but  no  indication  of 
a  separation  bubble  was  formd,  at  the  leading  edge.  The  prediction  for  the 
pressure  side  of  the  blade  matched  up  well  with  the  experimental  data.  For  the 
suction  side  of  the  blade,  the  shape  of  the  predicted  Cp  profile  seemed  to  agree  in 
the  axial  direction  but  the  magnitude  was  lower.  The  most  noticeable  difference 
was  that  the  diffusion  rate  for  the  CFD  analysis  from  approximately  0.45  to  0.7 
^/c  was  less  than  was  measured.  This  could  explain  why  there  was  no  boimdary 
layer  separation  predicted  at  the  midspan,  which  was  the  symmetry  plane  of  the 
computational  grid. 


Figure  27.  Predicted  vs.  Experimental  Cp  Distribution, 
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VI.  CONCLUSION  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

Compressor  stator  67B  cascade  blading  was  successfully  tested  at  an  off- 
design  inlet  air  angle  of  39.5°  in  the  Low-Speed  Cascade  Wind  Tunnel.  The 
experiments  were  conducted  at  a  inlet  test  section  Mach  number  of  .22  and 
Reynolds  number  of  640,000.  In  each  case,  total  plenum  pressure,  total  plenum 
temperature,  and  ambient  pressure  were  recorded  to  non-dimensionalize  all 
data.  The  tunnel  was  successfully  modified  by  inserting  a  fence  into  the  south 
endwall  boundary  layer  which  made  the  flow  more  symmetrical  about  the 
midspan  section.  This  made  the  midspan  LDV  measurements,  of  two  velocity 
components  in  coincidence  mode,  more  valid  than  previous  studies. 

Blade  surface  flow  visualization  was  successfully  performed  using  a 
titanium  dioxide  and  kerosene  mixture,  and  showed  the  symmetry  of  the  flow  at 
midpan  of  the  blades.  Periodicity  was  also  observed  with  this  technique. 

Blade  surface  pressure  measurements  were  obtained  and  compared  to 
both  previous  experiments  at  lower  inlet  flow  angles  and  numerical  predictions 
of  the  three-dimensional  flow  through  the  blade  row.  A  Cn  vs  p  curve  was 
generated  using  previous  data  from  on-design  incidence  at  pi=36.3°,  and  off- 
design  incidence  at  pl=38°,  and  data  from  the  current  study.  Results  in  Figure  28 
showed  that  the  blade  was  still  working  pi=39.5°,  even  though  boundary  layer 
separation  had  occurred. 
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34  36  38  40  42 

Beta 


Figure  28.  Cn  vs.  Beta  Curve. 

LDV  surveys  were  completed  at  8  different  stations  which  characterized 
the  flow  in  the  inlet,  the  blade  passage,  and  the  wake  area.  Reverse  flow  was 
measured  in  the  suction  surface  boundary  layer  aft  of  75%  axial  chord.  Reverse 
flow  was  also  measured  in  the  wake  at  105  and  120%  axial  chord  from  the 
leading  edge. 

CFD  analysis  was  performed  and  results  were  in  reasonable  agreement 
with  experimental  data.  Vortex  flow  at  the  trailing  edge  of  the  blade  and  near  the 
endwaU  was  indicated  by  the  solution.  However,  no  indication  of  reverse  flow 
was  found  at  midspan. 

B.  RECOMMENDATIONS 

Further  LDV  studies  should  be  performed  at  and  near  the  leading  edge  of 
the  blade  on  the  suction  side  to  see  if  a  separation  bubble  exists.  Three- 
dimensional  surveys  should  be  performed  in  order  to  characterize  the  flow  in 
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endwall  region.  This  will  allow  better  mapping  of  the  air  flow  to  more 
thoroughly  validate  the  CFD  code. 

Further  CFD  studies  should  be  initiated  to  try  and  match  the  Cp 
distribution  as  close  as  possible  by  varying  the  input  parameters.  Then,  and  only 
then  can  the  CFD  code  be  fully  validated  by  the  experimental  data. 
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APPENDIX  A. 


MIXING  INSTUCTIONS  FOR 
TITANIUM  DIOXIDE  &  KEROSENE 


STEPS 


1.  Put  6  Oz's  of  Pure  Vegetable  Oil  in  a  one-quart  plastic  container. 

2.  Put  about  1  tablespoon  of  Satiurn  yellow  pigment  in  sifter  and  sift  into 
container. 

3.  Mix  well  with  hand  stirrer/tongue  depressor. 

4-  Put  2  scoops  (little  plastic  cup)  of  TiOa  into  sifter  and  sift  as  much  as  you 
can  into  the  container.  Dump  the  big  chunks  back  into  Ti02  can. 

Repeat  this  step  2  more  times. 

Note:  It  may  take  10-15  minutes  for  each  time  of  sifting.  Mix  the  contents  by 
hand  before  each  repeat.  Do  not  attempt  to  push  the  Ti02  through  the 
sifter. 

5.  Add  4  Oz's  of  Kerosene  (more  can  be  added  to  make  mixtiue  thinner) 

6.  Add  20  squirts  of  oil  (SAE  30)  from  an  oil  can. 

7.  Put  on  magnetic  mixer  and  mix  for  10-15  minutes. 
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APPENDIX  B 


TABLE  OF  SCANIVALVE  PORTS 
AND 

CHANNEL  ASSIGNMENTS 


Scanivlave  #1 

Blade  Pressure  Measurements 


1 

25 

3  Suet.  Side 

2 

Calibration 

26 

4  Suet.  Side 

3 

Plenum  Press 

27 

5  Suet.  Side 

4 

18  Press  Side 

28 

6  Suet.  Side 

5 

17  Press  Side 

29 

7  Suet.  Side 

16  Press  Side 

30 

8  Suet.  Side 

Bi 

15  Press  Side 

31 

9  Suet.  Side 

8 

14  Press  Side 

32 

10  Suet.  Side 

9 

13  Press  Side 

33 

11  Suet.  Side 

10 

12  Press  Side 

34 

12  Suet.  Side 

11 

11  Press  Side 

35 

13  Suet.  Side 

12 

10  Press  Side 

36 

14  Suet.  Side 

37 

15  Suet.  Side 

38 

16  Suet.  Side 

7  Press  Side 

39 

17  Suet.  Side 

40 

18  Suet.  Side 

17 

5  Press  Side 

41 

19  Suet.  Side 

18 

4  Press  Side 

42 

43 

TE 

J0_ 

2  Press  Side 

44 

Blade  8, 1  Suet. 

45 

Blade  8, 2  Suet. 

LE 

46 

Blade  8, 3  Suet. 

23 

1  Suet.  Side 

47 

Blade  8, 4  Suet. 

24 

2  Suet.  Side 

48 

Blade  8, 5  Suet. 

Scanivalve  #2 

Rake  Probe  Measurements 


n 

25 

2 

Calibration 

26 

Rake  total 

3 

Plenum  Press 

27 

Rake  total 

4 

P  Wall  Statie 

28 

Rake  total 

5 

Not  Used 

29 

Rake  total 

Not  Used 

30 

Rake  total 

EH 

Not  Used 

31 

Rake  total 

8 

Not  Used 

32 

Rake  total 

9 

Not  Used 

33 

Rake  total 

10 

P  Prandtl  tot 

34 

Rake  total 

11 

P  Prandtl  stat 

35 

Not  Used 

12 

36 

Not  Used 

13 

Calibration 

37 

Not  Used 

14 

38 

Not  Used 

15 

Rake  total 

39 

Not  Used 

16 

Rake  total 

40 

Not  Used 

17 

Rake  total 

41 

Not  Used 

18 

Rake  total 

42 

Not  Used 

19 

Rake  total 

43 

Not  Used 

20 

Rake  total 

44 

Not  Used 

21 

Rake  total 

45 

Not  Used 

22 

Rake  total 

46 

Not  Used 

23 

Rake  statie 

47 

Not  Used 

24 

48 

Not  Used 

TABLE  -  B1 


TABLE  -  B2 
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APPENDIX  C 


FIND  (2-D)  SOFTWARE  INPUTS 

<A>  Color  link:  off 

Traverse:  TSl  Model  9500  Processors:  2  Mode:  Coincidence. 
Date  File:  d:\**\(filename)  Data  sample  size:  IK  Data  Points 

<I>  VO  PORT  AND  PROCESSOR  TYPE  SELECTION 

<P>  PROCESSOR  SETTINGS 

<0>  OPTICS  CONFIGURATION 

<E>  EXPERIMENT  DOCUMENTATION  AND  INPUTS 

<H>  HARDWARE  DIAGNOSTICS 

<F>  DATA  HLE  MANAGEMENT 

<T>  AUTOMATIC  TRAVERSE  PARAMETERS 

<R>  REALTIME  HISTOGRAM 

<M>  RETURN  TO  MAIN  MENU 

<F2>  ACQUIRE  DATA  FOR  (#  OF)  RAW  DATA  FILES 
<F3>  STORE  PROGRAM  DOCUMENTATION 

This  section  shows  what  is  tmder  each  sub-menu  and  the  inputs 

<I>  I/O  PORT  AND  PROCESSOR  TYPE  SELECTION 

Color  link:  off 

Traverse:  TSI  Model  9500  Processors:  2  Mode:  Coincidence. 

Date  File:  d:V*\(f Rename)  Data  sample  size:  IK  Data  Points 

VO  Port  Selections  LDV  Processor  Type 


Traverse  Controller 

=  COM2 

First  Processor 

=  1990 

Sony  Position  Encoder 

=  COMl 

Second  Processor 

=  1990 

Printer  Port 

=  LPTl 

Third  Processor 

=  1990 

Processor  VO 

=  COMl 

Color  Link 

=  Off 

Master  Interface 

=  1998A 

Program  Installation  Settings 

Computer  Bus  Type  =  PCBUS  Graphic  Type  =  EGA 

Monitor  Type  =  Color  Toggle  Selection  =  Highlight 

DMA  Chan:  1,  Port  Addr:  300H  Coml-Com4  Addr:  3F8H;2F8H;3E8H;2E8H 
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<P>  PROCESSOR  SETTINGS 


Number  of  Processors:  2 

Number  of  K-Data  Points:  1 K 

Data  Sampling  Method:  TBD-ON 

Coincidence  Window  width  (ps)  2.0  E5 

DMA  Timeout  (seconds)  300 

Acquisition  Mode  Coincidence 


Number  of  Cycles: 

Processor  Type 

Processor  Mode 

Filter  Range  High  limit: 

Low  limit: 

Timer  Comparison 
Gain 


Processor  1 
8 

1990 

CONT. 

50  Mhz 
.SMhz 
1 
1 


Processor  2 
8 

1990 

CONT. 

50  Mhz 
.3  Mhz 
1 
1 


<0>  OPTICS  CONFIGURATION 


Using  Half  Angle  Calculation 


Green  beam 

Blue  beam 

Fringe  Spacing  (microns) 

4.7569 

4.5119 

Frequency  Shift 

+5  Mhz 

+5  Mhz 

Half  Angle 

3.1 

3.1 

Focal  Length  (mm) 

762 

762 

Beam  Spacing  (mm) 

82.5 

82.5 

Wave  length  (nm) 

514.5 

488 
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APPENDIX  D 


LDV  SUMMARY  AND  REDUCED  DATA 


Station 

Survey 

Name 

Date 

Taken 

Re# 

xlOS 

Survey 

Points 

iH^gH 

mmm 

Tpi 

F 

1 

■iV4iUl=*M 

y05/97 

6.4 

84 

14.76 

12.1 

70 

77.034 

1 

EW.4I 

^07/97 

6.4 

84 

14.76 

12 

70 

76.724 

3 

0301inl3 

3/01/97 

6.4 

60 

14.78 

12 

70 

76.674 

6bl 

0211bl63 

2/11/97 

6.4 

40 

14.76 

11.9 

65 

76.051 

Yaw  4°  L 

7bl 

0211bl73 

2/11/97 

6.4 

40 

14.76 

11.9 

68 

76.268 

Yaw4'’L 

8bl 

6204bl8 

2/^97 

6.4 

47 

14.75 

12 

69 

76.677 

Yaw  4®  L 

8bl 

d207bl84 

2/97/97 

6.4 

47 

14.76 

12 

72 

76869 

Yaw4®L 

9bl 

0218bl93 

2/1S/97 

6.4 

60 

14.8 

12.2 

66.5 

76.370 

Yaw4®L 

llwk 

0302wkll 

l/m/97 

6.4 

n  50 

14.78 

12.3 

67.5 

77.417 

13wk 

0302wkl3 

3/02/97 

6.4 

52 

14.78 

12.2 

70 

77.292 
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Station  1  Inlet  survey 
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station  1 1nlet  Survey 
Vref  = 
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station  3  Inlet  survey 
Vref  = 

Blade  spacing  (s)  = 


0301  in  13 
76.6737  m/s 
152.4  mm 
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Station  6  Boundary  Layer  Survey 
Vref  = 

Blade  Chord  (c)  = 


0211bl63 
76.0511  m/s 
127.254  mm 
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station  7  Boundary  Layer  Survey  021 1  bl73 

Vref  =  76.2682  m/s 

Blade  Chorde  (c)  *  127.254  mm 
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station  8  Blade  #3  Boundary  Layer  Survey  0204bl83 

Vref  =  76.677  m/s 

Blade  Chord  (c)  =  127.254  mm 
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station  8  Blade  #4  Boundary  Layer  Survey 
Vref  = 

Blade  Chord  (c)  = 


0207bl84 
76.8688  m/s 
127.254  mm 


— .^..T.,..-«.t.-.T.T.._r,,,.,..  . 
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station  11  Wake  Survey  0302WK11 
Vref*  77.4166  m/s 

Blade  spacing  (s)  *  152.4  mm 
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station  13  Wake  Survey 
Vref  = 

Blade  spacing  (s)  = 


0302WK13 
77.2924  m/s 
152.4  mm 
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APPENDIX  E. 


STACK  AND  RVC3D  CODE  INPUTS 


Input  to  Stack 

&nll  km=70  rhub=0.0  rtip=0.998  nblade=l  ysp=0.0071  dhl=0.01 
dtl=0.30  &end 

Input  to  RVC3D 

'GELDER  CONTROLLED-DIFFUSION  CASCADE' 

&nll  im=340  jm=49  km=70  itl=80  iil=143  &end 

&nl2  cfl=5.0  aviscl=0.0  avisc2=0.0  avisc4=1.0  ivdt=l  nstg=4  itmax=7000 

irs=l  epi=.60  epj=.70  epk=.70  &end 

&nl3  ibcin=l  ibcex=l  isymt=l  ires=10  icmt=50 

iresti=0  iresto=l  ibcpw=0  iqin=0  &end 

&nl4  emxx=0.16976  emty=0.13994  emrz=0.0  expt=0.0  prat=0.9729  ga=1.4 

om=0.000000  igeom=0  alex=  9.0  &end 

&nl5  ilt=2  tw=1.00  renr=6.0e5  prm'=.7  prtr=.9  vispwr=.666666 

srtip=0.0  cmutm=10.  jedge=30  kedge=50  iltin=2  dblh=2.5  dblt=0.00  &end 

&nl6  iol=l  io2=340  oar=0  ixjb=0  njo=0  nko=0 

jo=0  5  ko=0  &end 
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APPENDIX  F. 


OUTPUT  FOR  INLET  AND  EXIT  CONDITIONS 


♦  J|c  **  Jjc  Jfj  He  *  sjc  ♦♦♦  Jlc  ♦  *************  He  ♦*♦****♦♦**  :|e  Jjc  ******  Jje  3|C  j|e  sfe  *  Jjc  Jfj  :|c  *  :Jc  3je  sjc  3|£  5|c  5|e  sjc  *  J|c  * 

j-direction  averaged  quantities  on  inlet 
derived  variables,  absolute  system 

rho  &  ps  are  area-averaged,  u,  v,  &  w  are  momentum-averaged  =>  approximate  mixed-out  average 
notation:  ii=rhoOref,  cr=cOref,  er=rr*cr**2,  pr=pOref,  tr=tOref,  alpha=atan(v/u),  phi=atan(w/u) 


k 

distance 

%  mdot 

vtot/cr  alpha 

phi 

ps/pr 

pO/pr 

ts/tr 

tO/tr 

Mach 

pO  loss 

1 

0.00000 

0.00000 

0.00000  0.00000 

0.00000 

0.96846 

0.96846 

1.00000 

1,00000 

0.00000 

0.00000 

2 

0.00267 

0.00034 

0.05901  56.02360 

0.00000 

0.96846 

0.97083 

0.99929 

0.99999 

0.05903 

0,00000 

3 

0.00541 

0.00128 

0.08180  46.49279 

0.00000 

0.96816 

0.97271 

0.99865 

0.99999 

0.08185 

0.00000 

4 

0.00821 

0.00265 

0.09594  41.84780 

0.00000 

0.96740 

0-97366 

0.99814 

0.99998 

0.09603 

0.00000 

5 

0.01109 

0.00432 

0.10526  39.79194 

0.00000 

0.96686 

0.97439 

0.99776 

0.99998 

0.10538 

0.00000 

6 

0.01403 

0.00619 

0.11125  38.99600 

0.00000 

0.96657 

0.97499 

0.99750 

0.99997 

0.11139 

0.00000 

7 

0.01706 

0.00821 

0.11513  38.83252 

0.00000 

0.96648 

0.97550 

0.99732 

0.99997 

0.11529 

0.00000 

8 

0.02016 

0.01034 

0.11778  38.96325 

0.00000 

0.96652 

0.97596 

0,99719 

0.99997 

0.11795 

0.00000 

9 

0.02334 

0.01257 

0.11979  39.19616 

0.00000 

0.96660 

0.97638 

0.99710 

0.99997 

0.11997 

0.00000 

10 

0.02660 

0.01488 

0.12152  39.42608 

0.00000 

0.96670 

0.97676 

0.99701 

0.99997 

0.12170 

0.00000 

11 

0.02996 

0.01728 

0.12315  39.60389 

0.00000 

0.96678 

0,97712 

0.99694 

0.99997 

0.12334 

0.00000 

12 

0.03340 

0.01977 

0.12477  39.71604 

0.00000 

0.96684 

0.97745 

0.99686 

0.99997 

0.12496 

0.00000 

13 

0.03694 

0.02237 

0.12639  39.76921 

0.00000 

0.96687 

0.97777 

0.99678 

0.99997 

0.12660 

0.00000 

14 

0.04058 

0.02506 

0.12802  39.77890 

0.00000 

0.96689 

0.97806 

0.99670 

0.99997 

0.12823 

0.00000 

15 

0.04432 

0.02787 

0.12961  39.76194 

0.00000 

0.96689 

0.97834 

0.99662 

0.99998 

0.12983 

0.00000 

16 

0.04817 

0.03080 

0.13115  39.73235 

0.00000 

0.96688 

0.97861 

0.99654 

0.99998 

0.13138 

0.00000 

17 

0.05213 

0.03385 

0.13263  39.69987 

0.00000 

0.96686 

0.97887 

0.99646 

0.99998 

0.13287 

0.00000 

18 

0.05621 

0.03703 

0.13404  39.67004 

0.00000 

0.96685 

0,97911 

0.99638 

0,99998 

0.13428 

0.00000 

19 

0.06042 

0.04033 

0.13538  39.64513 

0.00000 

0.96684 

0.97935 

0.99631 

0.99998 

0.13563 

0.00000 

20 

0.06475 

0.04378 

0.13665  39.62524 

0.00000 

0.96683 

0.97958 

0.99624 

0.99998 

0.13691 

0.00000 

21 

0.06921 

0.04736 

0.13788  39.60940 

0.00000 

0.96683 

0.97980 

0.99618 

0.99998 

0.13814 

0.00000 

22 

0.07382 

0.05109 

0.13905  39.59626 

0.00000 

0.96682 

0.98002 

0.99611 

0.99998 

0.13932 

0.00000 

23 

0.07857 

0.05497 

0.14019  39.58458 

0.00000 

0.96681 

0.98023 

0.99605 

0.99998 

0.14047 

0.00000 

24 

0.08348 

0.05901 

0.14129  39.57346 

0.00000 

0.96681 

0.98044 

0.99598 

0.99998 

0.14158 

0.00000 

25 

0.08856 

0.06321 

0.14237  39.56239 

0.00000 

0.96681 

0.98065 

0.99592 

0.99998 

0.14266 

0.00000 

26 

0.09380 

0.06760 

0.14342  39.55114 

0.00000 

0.96680 

0.98085 

0.99586 

0.99998 

0.14372 

0.00000 

27 

0.09922 

0.07216 

0.14445  39.53969 

0.00000 

0.96679 

0.98105 

0.99581 

0.99998 

0.14476 

0.00000 

28 

0.10483 

0.07692 

0.14546  39.52814 

0.00000 

0.96679 

0.98125 

0.99575 

0.99998 

0.14577 

0.00000 

29 

0.11064 

0.08188 

0.14646  39.51661 

0.00000 

0.96678 

0.98144 

0.99569 

0.99998 

0.14677 

0.00000 

30 

0.11666 

0.08706 

0.14744  39.50519 

0.00000 

0.96678 

0.98163 

0.99563 

0.99998 

0,14776 

0.00000 

31 

0.12290 

0.09247 

0.14840  39.49396 

0.00000 

0.96677 

0.98182 

0.99557 

0.99998 

0.14873 

0.00000 

32 

0.12937 

0.09811 

0.14935  39.48293 

0.00000 

0.96676 

0.98201 

0.99552 

0.99998 

0.14968 

0.00000 

33 

0.13610 

0.10401 

0.15029  39.47211 

0.00000 

0.96676 

0.98220 

0.99546 

0.99998 

0.15063 

0.00000 

34 

0.14308 

0.11018 

0.15122  39.46146 

0.00000 

0.96675 

0.98238 

0.99541 

0.99998 

0.15157 

0.00000 

35 

0.15034 

0.11663 

0.15214  39.45096 

0.00000 

0.96674 

0.98257 

0.99535 

0.99998 

0.15249 

0.00000 

36 

0.15789 

0.12339 

0.15305  39.44056 

0.00000 

0.96674 

0.98276 

0.99530 

0.99998 

0.15341 

0.00000 

37 

0.16576 

0.13047 

0.15396  39.43025 

0.00000 

0.96673 

0.98294 

0.99524 

0.99998 

0.15433 

0.00000 

38 

0.17396 

0.13789 

0.15487  39.42000 

0.00000 

0.96672 

0.98313 

0.99518 

0.99998 

0.15524 

0.00000 

39 

0.18251 

0.14568 

0.15577  39.40977 

0.00000 

0.96671 

0.98331 

0.99513 

0.99998 

0.15615 

0.00000 

40 

0.19144 

0.15386 

0.15667  39.39956 

0.00000 

0.96671 

0.98350 

0,99507 

0.99998 

0.15705 

0.00000 

41 

0.20076 

0.16246 

0.15757  39.38934 

0.00000 

0.96670 

0.98369 

0.99502 

0.99998 

0.15796 

0.00000 
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42  0.21052  0.17151  0.15847  39.37910  0.00000  0.96669  0.98388  0.99496  0.99998  0.15887  0.00000 

43  0.22074  0.18104  0.15937  39.36880  0.00000  0.96668  0.98407  0.99490  0.99998  0.15978  0.00000 

44  0.23146  0.19110  0.16027  39.35844  0.00000  0.96667  0.98426  0.99484  0.99998  0.16069  0^00000 

45  0.24270  0.20171  0.16118  39.34799  0.00000  0.96666  0.98445  0.99478  0.99998  0.16161  0.00000 

46  0.25451  0.21293  0.16210  39.33742  0.00000  0.96666  0.98465  0.99473  0.99998  0.16253  0.00000 

47  0.26694  0.22480  0.16303  39.32671  0.00000  0.96665  0.98485  0.99467  0.99998  0.16346  0.00000 

48  0.28004  0.23738  0.16396  39.31585  0.00000  0.96664  0.98505  0.99460  0.99998  0.16440  0.00000 

49  0.29385  0.25073  0.16490  39.30480  0.00000  0.96663  0.98525  0.99454  0.99998  0.16535  0.00000 

50  0.30844  0.26492  0.16586  39.29354  0.00000  0.96662  0.98546  0.99448  0.99998  0.16632  0.00000 

51  0.32388  0.28003  0.16683  39.28205  0.00000  0.96660  0.98567  0.99441  0.99998  0.16730  0.00000 

52  0.34025  0.29614  0.16782  39.27031  0.00000  0.96659  0.98589  0.99435  0.99998  0.16830  0.00000 

53  0.35763  0.31335  0.16883  39.25830  0.00000  0.96658  0.98611  0.99428  0.99998  0.16931  0.00000 

54  0.37611  0.33177  0.16985  39.24599  0.00000  0.96657  0.98634  0.99421  0.99998  0.17035  0.00000 

55  0.39580  0.35153  0.17091  39.23337  0.00000  0.96655  0.98658  0.99414  0.99998  0.17141  0.00000 

56  0.41684  0.37278  0.17199  39.22041  0.00000  0.96654  0.98682  0.99406  0.99998  0.17250  0.00000 

57  0.43936  0.39567  0.17309  39.20712  0.00000  0.96652  0.98707  0.99399  0.99998  0.17362  0.00000 

58  0.46352  0.42041  0.17424  39.19349  0.00000  0.96651  0.98733  0.99391  0.99998  0.17477  0.00000 

59  0.48951  0.44721  0.17542  39.17957  0.00000  0.96649  0.98760  0.99383  0.99998  0.17597  0.00000 

60  0.51754  0.47633  0.17665  39.16539  0.00000  0.96647  0.98789  0.99374  0.99998  0.17720  0.00000 

61  0.54787  0.50807  0.17792  39.15096  0.00000  0.96646  0.98818  0.99365  0.99998  0.17849  0.00000 

62  0.58080  0.54279  0.17925  39.13609  0.00000  0.96644  0.98849  0.99355  0.99998  0.17983  0.00000 

63  0.61666  0.58093  0.18065  39.12027  0.00000  0.96642  0.98882  0.99345  0.99998  0.18125  0.00000 

64  0.65587  0.62298  0.18213  39.10266  0.00000  0.96639  0.98917  0.99335  0.99998  0.18274  0.00000 

65  0.69892  0.66957  0.18370  39.08276  0.00000  0.96637  0.98955  0.99323  0.99998  0.18432  0,00000 

66  0.74641  0.72146  0.18536  39.06297  0.00000  0.96634  0.98995  0.99311  0.99998  0.18601  0.00000 

67  0.79906  0.77955  0.18709  39.05466  0.00000  0.96632  0.99038  0.99298  0.99998  0.18775  0.00000 

68  0.85776  0.84494  0.18879  39.08606  0.00000  0.96634  0.99085  0.99285  0.99998  0.18947  0.00000 

69  0.92361  0.91891  0.19030  39.20073  0.00000  0.96645  0.99136  0.99273  0.99998  0  19099  0  00000 

70  0.99800  1.00000  0.19080  39.23853  0.00000  0.96649  0.99153  0.99269  0.99998  0.19150  0.00000 

*****  overall  averages  ***♦*♦♦***♦********♦*****♦*♦**♦*****♦♦♦*♦♦**♦*♦♦♦♦♦**♦♦♦*♦♦*♦*♦*♦*♦* 

71  0.49900  0.15339  0.17298  39.22453  0.00000  0.96654  0.98706  0.99408  1.00007  0.17349  0.00000 

**************************************tHf*^;t************m***Mt*tf:1f*)Hf^:tl:*****m******** 

j-direction  averaged  quantities  on  exit 
derived  variables,  absolute  system 

rho  &  jjs  are  area-averaged,  u,  v,  &  w  are  momentum-averaged  =>  approximate  mixed-out  average 
notation:  n=rhoOref,  ci“c0ref,  ep=rr*cr**2,  pi=p0ref,  ti=t0ref,  alpha=atan(v/u),  phi=atan(w/u) 

k  distance  %mdot  vtot/cr  alpha  phi  ps/pr  pO/pr  ts/tr  tO/tr  Mach  pOloss 

1  0.00000  0.00000  0.00000  0.00000  0.00000  0.97290  0.97290  0.99999  0.99999  0.00000  0.02710 

2  0.00267  0.00031  0.03056  -7.45493  0.03446  0.97290  0.97354  0.99979  0.99997  0.03056  0.02646 

3  0.00541  0.00118  0.05509  -7.36759  0.00552  0.97290  0.97497  0.99937  0.99997  0.05511  0.02503 

4  0.00821  0.00252  0.07399  -7.25528  -0.01744  0.97290  0.97664  0.99889  0.99998  0.07403  0.02336 

5  0.01109  0.00425  0.08813  -7.11798  -0.03393  0.97290  0.97821  0.99845  1.00000  0.08820  0.02179 

6  0.01403  0.00629  0.09855  -6.96094  -0.05049  0.97290  0.97954  0.99808  1.00002  0.09864  0.02046 

7  0.01706  0.00858  0.10613  -6.78754  -0.06780  0.97290  0.98061  0.99779  1.00004  0.10624  0.01939 

8  0.02016  0.01109  0.11154  -6.59779  -0.08591  0.97290  0.98142  0.99756  1.00005  0.11168  0,01858 

9  0.02334  0.01376  0.11530  -6.38985  -0.10461  0.97290  0.98201  0.99740  1.00006  0.11545  0.01799 

10  0.02660  0.01658  0.11778  -6.16158  -0.12360  0.97290  0.98241  0.99729  1.00007  0.11794  0.01759 

11  0.02996  0.01954  0.11926  -5.91188  -0.14263  0.97290  0.98265  0.99722  1.00007  0.11943  0.01735 

12  0.03340  0.02260  0.11998  -5.64131  -0.16157  0.97290  0.98277  0.99718  1.00006  0.12015  0.01723 

13  0.03694  0.02577  0.12015  -5.35223  -0.18041  0.97290  0.98279  0.99717  1.00006  0.12032  0.01721 

14  0.04058  0.02902  0.11994  -5.04806  -0.19923  0.97290  0.98276  0.99717  1.00005  0.12011  0.01724 

15  0.04432  0.03237  0.11949  -4.73239  -0.21816  0.97290  0.98269  0.99718  1.00004  0.11966  0.01731 

16  0.04817  0.03581  0.11895  -4.40837  -0.23737  0.97290  0.98260  0.99720  1.00003  0.11912  0.01740 

17  0.05213  0.03934  0.11840  -4.07855  -0.25703  0.97290  0.98251  0.99722  1.00002  0.11857  0.01749 

18  0.05621  0.04296  0.11792  -3.74494  -0.27729  0.97290  0,98243  0.99723  1.00001  0.11808  0.01757 
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19  0.06042  0.04669  0.11754  -3.40912  -0.29828  0.97290  0.98237  0.99724  1.00001  0.11770  0.01763 

20  0.06475  0.05052  0.11729  -3.07238  -0.32011  0.97290  0.98233  0.99725  1.00000  0.11745  0.01767 

21  0.06921  0.05446  0.11717  -2.73562  -0.34284  0.97290  0.98231  0.99725  1.00000  0.11733  0.01769 

22  0.07382  0.05853  0.11716  -2.39916  -0.36648  0.97290  0.98231  0.99725  1.00000  0.11733  0.01769 

23  0.07857  0.06274  0.11726  -2.06255  -0.39101  0.97290  0.98232  0.99725  1.00000  0.11743  0.01768 

24  0.08348  0.06708  0.11745  -1.72450  -0.41638  0.97290  0.98235  0.99724  1.00000  0.11761  0  01765 

25  0.08856  0.07158  0.11769  -1.38302  -0.44252  0.97290  0.98239  0.99723  1.00000  0.11785  0  01761 

26  0.09380  0.07624  0.11797  -1.03604  -0.46937  0.97290  0.98244  0.99721  1.00000  0.11814  0.01756 

27  0.09922  0.08107  0.11828  -0.68187  -0.49685  0.97290  0.98249  0.99720  1.00000  0.11845  0.01751 

28  0.10483  0.08608  0.11860  -0.31928  -0.52490  0.97290  0.98254  0.99719  1.00000  0.11877  0.01746 

29  0.11064  0.09128  0.11893  0.05262  -0.55346  0.97290  0.98260  0.99717  1.00000  0.11910  0  01740 

30  0.11666  0.09668  0.11927  0.43440  -0.58245  0.97290  0.98265  0.99716  1.00000  0.11944  0.01735 

31  0.12290  0.10230  0.11960  0.82641  -0.61182  0.97290  0.98270  0.99714  1.00000  011977  0  01730 

32  0.12937  0.10813  0.11994  1.22879  -0.64149  0.97290  0.98276  0.99713  1.00001  0.12011  o'oi724 

33  0.13610  0.11421  0.12028  1.64157  -0.67136  0.97290  0.98282  0.99712  1.00001  0.12046  0.01718 

34  0.14308  0.12053  0.12064  2.06481  -0.70132  0.97290  0.98288  0.99710  1.00001  012082  0  01712 

35  0.15034  0.12712  0.12101  2.49835  -0.73126  0.97290  0.98294  0.99709  1.00002  0.12119  0.01706 

36  0.15789  0.13399  0.12140  2.94176  -0.76105  0.97290  0.98300  0.99707  1.00002  0.12158  0  01700 

37  0.16576  0.14116  0.12182  3.39423  -0.79060  0.97290  0.98307  0.99706  1.00002  0.12200  0.01693 

38  0.17396  0.14865  0.12226  3.85473  -0.81983  0.97290  0.98315  0.99704  1.00003  0  12244  0  01685 

39  0.18251  0.15648  0.12274  4.32197  -0.84869  0.97290  0.98323  0.99702  1.00003  0.12293  0  01677 

40  0.19144  0.16468  0.12326  4.79438  -0.87719  0.97290  0.98332  0.99699  1.00003  0.12345  0.01668 

41  0.20076  0.17327  0.12382  5.27007  -0.90535  0.97290  0.98341  0.99697  1.00004  0.12401  0.01659 

42  0.21052  0.18229  0.12442  5.74678  -0.93330  0.97290  0.98352  0.99694  1.00004  0.12461  0.01648 

43  0.22074  0.19177  0.12506  6.22180  -0.96119  0.97290  0.98363  0.99691  1.00004  0.12525  0.01637 

44  0.23146  0.20174  0.12574  6.69196  -0.98923  0.97290  0.98374  0.99688  1.00004  0  12594  0  01626 

45  0.24270  0.21224  0.12646  7.15372  -1.01773  0.97290  0.98387  0.99685  1.00005  0.12666  0  01613 

46  0.25451  0.22333  0.12722  7.60327  -1.04700  0.97290  0.98400  0.99681  1.00005  0.12742  0  01600 

47  0.26694  0.23504  0.12799  8.03682  -1.07750  0.97290  0.98414  0.99677  1.00005  0.12820  0.01586 

48  0.28004  0.24744  0.12878  8.45103  -1.10986  0.97290  0.98428  0.99673  1.00005  0.12899  0.01572 

49  0.29385  0.26058  0.12956  8.84368  -1.14499  0.97290  0.98442  0.99669  1.00005  0  12977  0  01558 

50  0.30844  0.27453  0.13031  9.21451  -1.18379  0.97290  0.98455  0.99665  1.00005  0.13052  0^01545 

51  0.32388  0.28934  0.13100  9.56664  -1.22868  0.97290  0.98468  0.99662  1.00005  0.13123  0.01532 

52  0.34025  0.30509  0.13159  9.88626  -1.28193  0.97290  0.98478  0.99659  1.00005  0  13182  0  01522 

53  0.35763  0.32185  0.13201  10.14582  -1.34402  0.97290  0.98486  0.99657  1.00005  0.13224  0.01514 

54  0.37611  0.33967  0.13225  10.32509  -1.41326  0.97290.0.98490  0.99656  1.00006  0.13247  0.01510 

55  0.39580  0.35865  0.13233  10.41270  -1.48424  0.97290  0.98492  0.99657  1.00007  0.13256  0.01508 

56  0.41684  0.37887  0.13236  10.40686  -1.54783  0.97290  0.98492  0.99658  1.00008  0.13258  0  01508 

57  0.43936  0.40046  0.13242  10.31474  -1.59170  0.97290  0.98494  0.99659  1.00009  0.13265  0.01506 

58  0.46352  0.42360  0.13265  10.15019  -1.60214  0.97290  0.98498  0.99659  1.00011  0.13288  0.01502 

59  0.48951  0.44851  0.13315  9.93015  -1.56696  0.97290  0.98507  0.99657  1.00012  0.13338  0.01493 

60  0.51754  0.47552  0.13398  9.67118  -1.47887  0.97290  0.98522  0.99653  1.00012  0.13421  0  01478 

61  0.54787  0.50504  0.13519  9.38806  -1.33798  0.97290  0.98545  0.99646  1.00012  0.13543  0.01455 

62  0.58080  0.53759  0.13682  9.09526  -1.15267  0.97290  0.98576  0.99636  1.00010  0.13707  0  01424 

63  0.61666  0.57381  0.13891  8.80923  -0.93853  0.97290  0.98615  0.99622  1.00008  0.13917  0.01385 

64  0.65587  0.61444  0.14147  8.54851  -0.71573  0.97290  0.98665  0.99606  1.00006  0.14175  0.01335 

65  0.69892  0.66030  0.14442  8.33010  -0.50528  0.97290  0.98724  0.99587  1.00004  0.14472  0  01276 

66  0.74641  0.71231  0.14757  8.16412  -0.32497  0.97290  0.98788  0.99567  1.00003  0.14789  0.01212 

67  0.79906  0.77148  0.15060  8.05064  -0.18596  0.97290  0.98850  0.99548  1.00002  0.15094  0.01150 

68  0.85776  0.83886  0.15312  7.98126  -0.09066  0.97290  0.98904  0.99533  1.00002  0.15348  0  01096 

69  0.92361  0.91561  0.15479  7.94465  -0.03275  0.97290  0.98940  0.99523  1.00002  0.15516  0  01060 

70  0.99800  1.00000  0.15536  7.93298  0.00000  0.97290  0.98952  0.99519  1.00002  0.15574  0.01048 

*****  overall  averages  *********************♦***♦**♦***♦♦***♦♦**♦♦♦*♦♦*♦♦♦*♦♦♦*♦♦♦♦*♦♦****♦*** 

71  0.49900  0.15346  0.13682  7.28813  -0.69505  0.97290  0.98576  0.99634  1.00008  0.13707  0.01424 
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